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TRANSIENT  HEATING  OF  THIN  PUTBS 

ABSTRACT 

Two  classes  of  heating  problems  involving  the  heating  of  thin  plates 
over  a  portion  of  their  surface  are  considered.  Such  problems  have 
mumerous  applications,  one  such  application  being  the  determination  of 
the  temperature  rise  in  materials  upon  which  a  laser  beam  is  directed. 

The  heating  effects  caused  by  a  laser  beam  fits  into  the  class  of 
problems  for  which  the  heating  source  has  circular  symmetry.  The  second 
class  of  problems  arises  when  the  source  term  has  symmetry  about  a  line 
on  the  surface.  Solutions  in  the  form  of  a  definite  integral  are 
obtained  for  both  classes  of  problems  for  a  general  time  independent 
source  of  unspecified  spatial  distribution.  Several  solutions  in  closed 
form  are  obtained  for  particular  source  functions  and  for  particular 
values  of  the  various  parameters  of  the  problem.  For  the  case  of  a  disk 
source  an  error  bound  useful  in  the  numerical  quadrature  of  the  integral 
solution  is  obtained. 
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dimensionless  heat  transfer  coefficient.  See  Eq. (8) 
temperature 
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I.  INTRODUCTION 


Several  types  of  conduction  heat  transfer  applications  involve  the 

determination  of  transient  temperature  in  a  thin  plate  heated  over  a 

portion  of  its  surface,  the  surface  being  simultaneously  cooled  through 

surface  convection.  Typical  applications  include  the  heating  of  thin 

metal  plates  with  a  laser  or  incendiary,  drilling  holes  in  metal, 

acetylene  burniig,  arc  welding  and  friction  heating  of  machine  parts. 

Solutions  for  uniform  disk  and  unifonn  strip  sources  were  obtained  by 
1  ^ 

Thomas  ,  the  solutions  being  in  the  form  of  definite  integrals  for  the 
most  general  cases.  For  special  cases  such  as  no  cooling,  steady  state, 
or  at  the  center  of  the  heated  area,  Thomas  reduced  the  general  integral 
solutions  to  the  form  of  functions  for  which  tables  exist.  In  this 
report  the  results  cited  above  are  generalized  to  include  solutions  for 
all  time  independent  source  functions  meeting  certain  fairly  general 
criteria.  Additionally  several  new  results  in  closed  form  are  obtained. 
For  the  case  of  a  disk  source  the  general  integral  solution  is  in  the 
form  of  a  definite  integral  on  the  interval  (0,«).  For  this  case  an 
error  bound  is  obtained  which  yields  an  upper  bound  on  the  error  in¬ 
curred  when  this  integral  is  approximated  by  numerical  quadrature  over 
finite  limits. 


II.  CIRCULAR  SYMMETRIC  SOURCES 
We  assume  the  thin  plate  of  thickness  &  has  constant  thermal 
properties,  the  cooling  effect  can  be  represented  by  Newton's  Law  of 
Cooling  and  the  temperature  gradient  across  the  thickness  of  the  thin 
plate  is  negligible.  We  also  assume  that  the  initial  temperature  and 
the  temperature  of  the  cooling  medium  is  constant  at  =  0^.  Under 
these  conditions  and  for  a  circular  symmetric  source.  Fig.  1,  the 
applicable  partial  differential  equation,  initial  condition  and 
boundary  conditions  are  the  following; 

„  ^(0  -  0  )  +  R  >  0,  T  >  0 


pc 


R 


0 


o 

0. 


R  ^  0,  T  =  0 


(1) 

(2) 


Referenaea  are  Hated  on  page  24. 
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CIRCULAR  SYMMETRIC  SOURCE 


h  (e-do) 


Fig.  1 
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<  «  R 

In  the  usual  fashion  it  is  desirable  to  work  with  dimensionless 
variables.  If  we  let 


®,  T  is  0 

(3) 

0,  T  ^  0. 

(4) 

t  =  qt/3'^. 

(5) 

r  =  RiS'S 

(6) 

T  =  k6(fl  - 

(7) 

=  2h/3^k'^6"\ 

(8) 

q(r)  =  Q(j3r,fl)Q^^, 

(9) 

where 


%  -  Q(0,i3), 

1,^-1  -1 
a  =  kp  c 


(10) 


(11) 


then  the  problem  can  be  expressed  in  the  following  dimensionless  form; 


3T  rx  .  1  3T 


(12) 


at  -  +  •  «T  +  <)W>  r>0.  t  >0 

T  =  0,  r'^0,t  =  0  (13) 

r  -*  ®,  t  ^  0  (14) 


T  =  =  0, 


T<«  r  =  0,  t^O.  (15) 

The  technique  employed  to  obtain  a  solution  to  (12)  -  (15)  is  the 
integral  transform  method.  The  zero  order  Hankel  transform  of  any 
arbitrary  function  v(r)  is  defined  by 

Vj^(o)  =  J  rv(r)J^(ar)dr,  (16) 


where  its  inverse  transform  is  given  by 


00 

v(r)  =  J  avjj(a)J^(ar)da. 


(17) 


By  taking  the  zero  order  Hankel  transform  of  (12)  ,  we  obtain  the 
subsidiary  equation 


11 


(o I t)  n  n 

- =  -(a^  +  €^)Tj^(a,t)  +  q^(a).  (18) 

2 

This  result  arises  from  the  identity  (see  Carslaw  and  Jaeger  ) 

o  o 

which  holds  if  v(0)  is  finite  and 

v(r)  »  -♦  0  as  r  -*  ».  (20) 

or 

From  (13)  and  (16)  we  have 

Tj^(o,0)  -  0.  (21) 

The  solution  of  (18)  subject  to  (21)  is 

fj^(o,t)  =  qj^(o)(a^  +  €^)’^[1-  exp[-t(o^  +  f^)]].  (22) 

From  (22)  and  (17)  we  have 

oo 

T(r,t)  =  J  qjj(a)(a^  +  C^)’^0'JQ(C7r)Cl-exp[-t(a^+C^)  ]]da.  (23) 

o 

This  solution  is  valid  for  all  circular  S3rmmetric  sources  for  which  a 
Hankel  transform  q„(o)  is  defined.  For  the  most  general  cases  the  above 
integral  cannot  be  reduced  further  and  one  must  resort  to  numerical 
quadrature  for  obtaining  numerical  values. 

A.  Gaussian  Circular  Source 

A  Gaussian  circular  source  is  defined  by 

2  -2 

Q(R,i3)  -  ^  ,  (24) 

which  in  dimensionless  form  reduces  to 

2 

q(r)  »  e  .  (25) 

We  consider  the  case  for  no  cooling,  c  *  0.  For  this  case  ^„(C7)  is 

H 

★  3 

given  by  (10),  p.29,  Brdelyi  ,  et.  al  ,  Vol.2 

Hereafter  the  two  volumes  of  this  reference  will  be  referred  to  as 
[IMOTI]  and  [BM0T2]. 
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t 


(26) 


(27) 


qj^(o)  = 

The  general  solution  (23)  reduces  to 

T(r,t)  =  ^  J  a  ^J^(or)j^e  ^  -  e 
*^0 

To  evaluate  this  Integral  we  make  use  of  the  identity 

0  ^e  J^(ar)do  =  ^(ff/a)^  J  ®  ^  Ii[®^(8a)  ^]do, 

‘^o  ‘^r  * 

which  i.s  derived  by  combining  (7)  p.5  and  (8)  p.29  of  [BM0T2].  We  also 
have  the  identity 


(28) 


I^(z)  =  (2ir2)”^(e^  -  e  ^). 


Hence 


00  2  ® 

a  ^e  J  (ar)do  «  ra^(l-e^ 

J  O  J 


2  -1 
a 


)da. 


(29) 


(30) 


2  "1 

By  letting  y  =  ^  a  ,  wa  obtain  the  following  identity: 


*«  .  1  2  -1 


da-i 


-y  -1 

2  -1®  y 

a 


From  (27) ,  (30)  and  (31)  we  obtain  the  final  result 

T(r,t)  -  i{B^[r^(4t+l)"^]  -  B^(r^)}. 


(31) 


(32) 


We  note  that  this  solution  is  unbounded  for  r  »  0;  thus, we  treat  this 
case  seperately.  For  r  =  0  in  (27),  we  obtain 


-y  -1 


T(O.t)  "  C  J  2  e  'y  'dy  - 


00 


-y  -1 
»  «» 


,e  'y 


dy] 


2 

2 


-y  -1, 
e  dy. 


.  Um  i  r'^ 

‘  JrV(4t+l) 

Since  the  integral  is  to  be  evaluated  near  the  origin,  we  can  let 

e  ^  =  1  -  y  +  0(y^ 


(33) 


(34) 


13 


obtaining 


2 

T(0,t)  =.  J  2  [(1  -  y)y'^  +  0(y)]dy  (35) 

/  (4t+l) 

-  iln(4t+  1)  +^'t'”o[‘>(y'>lr'/(4t+l)] 

»  ^In  (4t  +  1) . 

Hence,  the  solution  for  a  Gaussian  circular  source  and  no  cooling  is 
given  by 


T(r,t) 


i(E2[r^4t+l)'^]  -  EjCr^), 
iln(4t  +  1), 


r  >  0 
r  »  0. 


(36) 


B.  Uniform  Disk  Source 

A  uniform  disk  source  is  defined  by 

Qq.  R  ^  ^ 

Q(R,fl)  - 

0,  R  >  jS, 

which  in  dimensionless  form  reduces  to 

U,  r  ^  1 
*jo,  r  >1. 

The  zero  order  Hankel  transform  of  (38)  is  given  by  (3)  ,  p.47  of 
[BMOT2]. 


(37) 


(38) 


qj^(a)  -  J^(C7)/a.  (39) 

For  this  case  the  general  solution  reduces  to 

GO 

T(r,t)  -  J  J^(ar)J2((r)<o^  +  -  exp[-t(o^  +  c^}]da.  (40) 

o 

an  integral  solution  obtained  by  Thomas^.  In  this  form  the  above 
integral  cannot  be  reduced  further.  For  steady  state  the  integral 
reduces  to 


T(r,«)  »  J  J^(ar)Jj^(o)(a^  + 

*^o 

For  r  >  1  this  integral  is  given  by  (12),  p.49 


C^'^da.  (41) 

of  [BM0T2],  namely. 
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T(r,®)  =  €"4j^(c)K^(€r)  ,  r  >  1.  (42) 

For  r  <  1  we  delete  the  transient  term  in  the  original  partial 
differential  equation  and  obtain 

+  1  -  0.  (43) 

The  homogeneous  equation  is  one  of  the  standard  Bessel  types  which  has 

solutions  I  (cr)  and  K  (cr).  We  can  eliminate  K  (cr) ,  since  we  expect  a 
o  o  o 

finite  temperature  at  r  =  0.  The  nonhomogeneous  equation  thus  has  a 
solution 

T(r,®)  =  f  ^  +  Cl^(cr),  r  <  1.  (44) 

We  determine  the  constant  C  by  forcing  the  solution  to  be  continuous  at 
r  =  1,  that  is,  equating  (42)  and  (44)  at  r  =1.  We  obtain 

c  >  (45) 

By  employing  tht  Wronskian 

lj(OK^(0  +  I„(e)K^(e)  -  (46) 

we  have 

C=-€"^K^(C).  (47) 

This  yields 

T(r,«)  =  (48) 

The  steady  state  solution  for  a  disk  source  with  cooling  is  thus 
given  by 

C  ^[1  -  CK,(Ol  (€r)],  r  i  1 

T(r,co)  »  °  (49) 

€  ^I,  (€)K  (cr)  ,  r  >  1. 

1  o 

For  the  general  case  (40)  can  be  approximated  by  an  integral  over 
finite  limits.  Let 

T(r,t)  =  J  J^(ar)J^(a)(o^  +  €^’^[1  -  exp{-t(a^  +  €^)}]do  +  B.  (50) 

'*0 
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We  will  now  determine  a  U  such  that  B  is  bounded  by  some  a  priori 
acceptable  error.  First  we  determine  a  bound  for  and  J^(z)t  From 

Watson^,  p.74 


J^(z)  »  +  H^^^(z)], 


(51) 


and  Watson  p.  168 


H^^^(z)  «  (2/ffz)^e^®/r(n+^)  f  e  u”’^(l  +  ^iu/z)””^du, 
n  J 


»  (2/irz)^e  ^'^/r(n+^)  '  e  ”  u""^(l  -  ^iu/z)"'^du,  (52) 


n 


where  Q  z  -  ntr/2  -  tr/4,  and  2n  is  not  an  odd  positive  integer.  Since 
|l  -  iu/22|  »  |l  +  iu/2z|  a  [l  +  (u/2z)^]^,  (53) 

then 

|j^(z)|^  (2/frz)^[r(n+^)  ]  ^  e  '^j^uCl  +  (u/2z)^]^J*'  du.  (54) 

For  n  »  0,  we  have 

|j^(z)|  s  (2/irz)^[r(i)]’^J  u'^e'“[l  +  (u/2z)^r^du 


s  (2/ire)^[ni) I'H  u^'*e'“du  -  (2/m)^.  (55) 

o 

For  n  a  I,  where  we  assume  z  we  obtain 

|jj(z)|  fi  (2/trz)^[r(i)]'^J  u^[l  +  (u/2z)^]S'“du 


S  (2/lFz)^[r(|)  ]  [u^  +  u*]*e  “du 


The  bounds  derived 
forth  in  more  detail 


here  are  generalized  and  the  derivation  is  put 
in  a  forthcoming  report  by  the  authors. 
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I 


^  (2/ffz)^[n|)]‘^[  J  2S’“ciu  +  J^a^e'^'clu] 


Now  in  (50) 


7  -1  -1 

^  2*(1  +  e  )ir  (z)T 


00 

|e|  s  f  |J  (ar)|  lj,(a)|  (o^  +  tS'^da. 

Jy'  O  1 


(56) 


(57) 


If  we  apply  the  bounds  (55)  and  (56)  to  (57) ,  this  yields  (assuming 
r  >  0) 


.03 

E|  ^  ABr*  [a(a^  +  c^)]  do. 


'U 


where 


and 


i  -1  -1 
A  =  2^(1  +  e  ^)ff  " 


B  -  (2/ff)^. 


(58) 


(59) 


(60) 


This  implies  that 

|e|  s  ^ABr'^t'^ln[(U^  +  (61) 

If  we  solve  for  U,  we  obtain  the  limit  of  integration,  that  is, 

U  -  cCexp{2rV|B|A*^B“^}  -  1]"^,  (62) 

where  E  is  the  error.  For  the  case  when  r  »  0,  we  use  the  fact 


J  (0)  »  1.  Hence 
o 


I  *4” 


■i«-2 


a  ®o  do 


i  fAU"^, 


or 


U  -  (§A|B|"S^ 


(63) 


(64) 


III.  SOURCES  SYMMETRICAL  ABOUT  A  LINE 
The  partial  differential  equation,  initial  and  boundary  conditions 
governing  this  case  in  dimensionless  form  are 


^  +  <,(x).  X  >0.  t  >0 

3x 


(65) 
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T  =  0, 

X  ^  0, 

t 

=  0 

(66) 

Bx 

X  ■  0, 

t 

>  0 

(67) 

O 

tl 

II 

H 

X  -►  00, 

t 

^  0. 

(68) 

Condition  (67)  arises  from  the  assumption  that  the  heat  source  is 
symmetrical  about  the  line  x  »  0,  that  is, 

q(x)  =  q(-x).  (69) 

Again  we  employ  the  integral  transform  method  for  obtaining  a  solution. 
We  define  the  Fourier  cosine  transform  of  an  arbitrary  function  u(x)  as 


u^(o)  =  (2/tr)^J  u(x)cos(ax)dx, 
o 

and  its  inverse  transform  by 

u(x)  »  (2/7r)^J  u^(cr)cos(ax)do. 
Employing  the  identity 

1  ®  2  T 

(2/ff)®J  ^-|cos(ax)dx  »  -a^u^(a)  -  (2/ff)®|^ 
odx 


(70) 


(71) 


XaO 


(72) 


and  taking  the  Fourier  cosine  transform  of  (65)  ,  we  obtain  the  subsidiary 
equation 


dT^(a,t) 

dt 


(0^  +  €^)f^(a,t)  +  q^(o), 


(73) 


which  is  subject  to 


f  (0,0)  -  0.  (74) 

c 

As  before  the  solution  is 

f^(o,t)  »  q^(a)[o^  +  c^]’^[l  -  exp{-t(a^  +  c^)}].  (75) 

Taking  the  inverse  Fourier  cosine  transform,  we  obtain 
-  00 

T(r,t)  -(2/Tr)®J  q^(o)[0^  +  €^]’^[l  -  exp{-t(a^  +  c^)  ]co8(xo)da.  (76) 

This  solution  is  valid  for  all  sources  S3mm]etrical  about  the  line  x  »  0, 
for  which  a  Fourier  cosine  transform  exists  and  for  which  (69)  holds. 
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A.  Uniform  Strip  Source 

A  uniform  strip  source  Is  defined  by 


Q(x,/3) 


jv 

(o,  X  >  /3. 


(77) 


In  dimensionless  form  this  equation  becomes 


q(x) 


(1,  X  SI 

10.  X  >  1. 


The  Fourier  cosine  transform  of  (78)  Is 

q  (o)  =  (2/ir)^a  ^sln(o). 


(78) 


(79) 


Inserting  (79)  Into  (76)  ,  we  obtain  the  Integral  solution  obtained  by 
Thomas^ 

T(x,t)  =  (2/77)J  a  ^sln(o)[o^  +  “  exp{-t(a"  +  C^)}> 

*^0 

xCOs(xo)dff.  (80) 


Let 


where 


T(x,t)  »  Tj^(x)  -  T2(x,t)  , 


(81) 


and 


00 

Tj^(x)  »  (2/ir)J  0  ^sln(a)[o^  +  e^]  ^  cos(xp)do  (82) 


2  00  2 

T2(x,t)  =  (2/ff)e  ^  <7  ^sln(a)[a^  +  C^]'^e  cos(xa)da.  (83) 


From  (4),  p.l9,  of  [EMOTI] 


i  €'"[1  - 

•l^*)  "j  -2  -€x 

Ice  t 


-  e  ^co8h(€x)],  0  ^  X  ^  1 


8lnh(c)  ,  1  <  X. 

In  evaluating  (83)  ,  we  first  consider  the  related  Integral 


(84) 


2  „oo 


w(x,t)  ■  (2/ir)e  ^  ^  a  +  €^]  ^ 

V 


e  sln(xa)da 


(85) 


(2/ir)e'^e’“^J  [a'^  -  9(9^  + 


[e  sln(xa)do.  (86) 
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By  employing  (21),  p,73  and  (26),  p.74  of  [BMOTl],  we  obtain 

2 

w(x,t)  »  €  erf(^xt"^)  + 

^[e^*erfc(ct^  +  jtxt”^)  -  e  ^*erfc(ct^  -  ^xt”^)].  (87) 

It  can  be  shown  that 

T2(x,t)  -  ^[w(l+x,t)  +  w(l-x,t)].  (88) 

If  we  collect  these  results,  the  solution  for  a  unifoinn  strip  source  is 
given  by 

T(x,t)  -  T^(x)  -  T2(x,t),  (89) 


where 


and 

T2(x,t) 


C  ^[1  -  e  ^cosh(cx)],  0  ^  X  1 

Ti(x)  =■  .2  -€X 

C  e  sinh(€) ,  1  <  x 

2 

[erf{i(l+x)t’^}  -  erf {^(l-x)t"^}J  + 


(90) 


(91) 


^^e^^^'*'*^erfc{€t^  +  ^(l+x)t'^}  -  e"^^^'*’*^erfc{ct^  -  ^(l+x)t”^}  + 
gC(l  *^erfc{€t^  +  ^(l-x)t"^}  -  e  *^erfc{et^  -  ^(l-x)t 
B.  GAUSSIAN  STRIP  SOURCE 


We  define  a  "Gaussian  strip  source"  by 

-x2/fl2 

<3(x./J>  =■  Q^e  *  . 

In  dimensionless  form  this  source  is  given  by 

_  2 

q(x)  -  e’* 

and  its  Fourier  cosine  transform,  which  is  obtained  from  (11), 
p.  15  of  [KMOTl]  is 


(92) 


(93) 


2-^e’" 


(94) 
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The  solution  for  a  strip  source  is  thus 

T(x,t)  =  e’^  [a^  +  €^l"^[l  -  exp{-t(ff^  +  €^) )  ]co8(ax)da.  (95) 

Employing  (15),  p.  15  of  [EMOTl]  we  obtain 

2 

T(x,t)  =  |e  ^^[erfc(^€-x)  - 

erfc{^€(4t+l)^  -  x(4t+l)"^}]  + 

e^^[erfc(^C+x)  -  erfc{^e(4t+l)^  +  x(4t+l)  "^]  ]j-.  (96) 

Note  that  in  the  solutions  for  sources  symmetrical  about  a  line,  the 
domain  of  interest,  -®  <  x  <  ®,  has  been  restricted  to  0  ^  x  <  ®  for 
mathematical  simplicity.  This  causes  no  loss  of  generality,  since  the 
solutions  obtained  can  be  extended  to  negative  x  simply  by  employing 
the  symmetry  condition  T(-x)  »  T(x),  x  ^  0, 

IV.  GRAPHICAL  RESULTS 

Figure  2  is  a  graphical  presentation  of  the  dimensionless 
temperature  rise  at  the  center  of  the  heated  uniform  strip  divided  by 
the  dimensionless  time  as  a  function  of  the  dimensionless  group  /3(aT) 
By  use  of  the  dimensionless  variables  defined  in  the  List  of  Sjnnbols 
the  effect  on  the  temperature  rise  above  ambient,  (6-0^),  can  be 
determined  as  a  function  of  any  of  the  variables  k,  6,  h,  Q^,  o',  T* 

Figure  3  is  a  graph  of  the  same  quantities  as  Figure  2  for  no 
surface  cooling  and  is  a  comparson  of  the  transient  temperature  rise 
above  ambient  at  the  source  center  for  the  sources  listed. 
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Fig.  3.  Effect  of  source  size,  source  type,  and  heating  time  on 
maximum  temperature  of  thin  plate  with  no  cooling  (C"0) . 
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Washington,  D.C, 

l>.  ARSTRACT 

Two  classes  of  heating  problems  involving  the  heating  of  thin  plates  over  a  portion 
of  their  surface  are  considered.  Such  problems  have  numerous  applications,  one  such 
application  being  the  determination  of  the  temperature  rise  in  materials  upon  which  a 
laser  beam  is  directed.  The  heating  effects  caused  by  a  laser  beam  fits  into  the  class 
of  problems  for  which  the  heating  source  has  circular  symmetry.  The  second  class  of 
problems  arises  when  the  source  term  has  symmetry  about  a  line  on  the  surface. 

Solutions  in  the  form  of  a  definite  integral  are  obtained  for  both  classes  of 
problems  for  a  general  time  independent  soxirce  of  unspecified  spatial  distribution. 
Several  solutions  in  closed  form  are  obtained  for  particular  source  functions  and  for 
particular  values  of  the  various  parameters  of  the  problem.  For  the  case  of  a  disk 
source  an  error  bound  useful  in  the  numerical  quadrature  of  the  Integral  solution  is 
obtained. 


